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We use observational data from Type la Supernovae (SNIa), Baryon Acoustic Oscillations (BAO), 
and Cosmic Microwave Background (CMB), along with requirements of Big Bang Nucleosynthesis 
(BBN), to constrain the cosmological scenarios governed by Hofava-Lifshitz gravity. We consider 
both the detailed and non-detailed balance versions of the gravitational sector, and we include the 
matter and radiation sectors. We conclude that the detailed-balance scenario cannot be ruled out 
from the observational point of view, however the corresponding likelihood contours impose tight 
constraints on the involved parameters. The scenario beyond detailed balance is compatible with 
observational data, and we present the corresponding stringent constraints and contour-plots of the 
parameters. Although this analysis indicates that Hofava-Lifshitz cosmology can be compatible 
with observations, it does not enlighten the discussion about its possible conceptual and theoretical 
problems. 
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I. INTRODUCTION 

Recently, a power-counting renormalizablc, ultra- 
violet (UV) complete theory of gravity was proposed by 
Hof ava in [JJ, d, H, H[ . Although presenting an infrared 
(IR) fixed point, namely General Relativity, in the UV 
the theory possesses a fixed point with an anisotropic, 
Lifshitz scaling between time and space. Due to these 
novel features, there has been a large amount of effort 
in examining and extending the properties of the the- 
ory itself EJl.JB H & {$L El M iiEi 
El EI EI S3, US, M M M M M M KM- Addi- 
tionally, application of Hofava-Lifshitz gravity as a cos- 
mological framework gives rise to Hofava-Lifshitz cos- 
mology, which proves to lead to interesting behavior 
[H |3Q |. In particular, one can examine specific solu- 
tion subclasses |Hj M, lajiJilM H E3, EH, M, El , 

J , the gravitational wave 
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l81|, the thermo- 
However, despite 



this extended research, there are still many ambiguities 
if Hofava-Lifshitz gravity is reliable and capable of a suc- 
cessful description of the gravitational background of our 
world, as well as of the cosmological behavior of the uni- 
verse fTTl [TH [TM [U " 



Although the discussion about the foundations and the 
possible conceptual and phenomenological problems of 
Hofava-Lifshitz gravity and cosmology is still open in 



'Electronic address: 
t Electronic address: 



sourish.d@gmail.com 



the literature, it is worth investigating in a systematic 
way the constraints imposed by observations in a uni- 
verse governed by Hof ava gravity. Thus, in the present 
work we use Big Bang Nucleosynthesis conditions, to- 
gether with Type la Supernovae (SNIa), Baryon Acoustic 
Oscillations (BAO) and Cosmic Microwave Background 
(CMB) data, in order to construct the corresponding- 
probability contour-plots for the parameters of the the- 
ory. Furthermore, in order to be general and model- 
independent, we perform our analysis with and without 
the detailed-balance condition. As we will show, both 
the detailed-balance and beyond-detailed-balance formu- 
lations are compatible with observations, however under 
tight constraints on the model parameters. 

The paper is organized as follows: In section QTJ we 
present the basic ingredients of Hofava-Lifshitz cosmol- 
ogy, extracting the Friedmann equations, and describ- 
ing the dark matter and dark energy dynamics. In sec- 
tion IIIII we constrain both the detailed-balance and the 
beyond-detailed-balance formulations from the observa- 
tional point of view. Finally, section HVl is devoted to the 
summary of the obtained results. 



II. HORAVA-LIFSHITZ COSMOLOGY 

Let us present the scenario where the cosmological evo- 
lution is governed by Hofava-Lifshitz gravity [H [30j ] . 
The dynamical variables are the lapse and shift func- 
tions, N and Ni respectively, and the spatial metric gij 
(roman letters indicate spatial indices) . In terms of these 
fields the full metric is written as: 

ds 2 = -N 2 dt 2 + g lJ {dx l + N l dt)(dx 3 + N j dt), (1) 

where indices are raised and lowered using . The scal- 
ing transformation of the coordinates reads: 
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A. Detailed Balance 

The gravitational action is decomposed into a kinetic 
and a potential part as S g = J dtd 3 x y /gN(C K +C v )- The 
assumption of detailed balance 0] reduces the possible 
terms in the Lagrangian, and it allows for a quantum 
inheritance principle [l|, since the (D + l)-dimensional 
theory acquires the renormalization properties of the D- 
dimensional one. Under the detailed balance condition 
the full action of Hofava-Lifshitz gravity is given by 



S g = 



dtd 3 Xy/gN I —{K^K 10 - XK 2 ) 
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is the extrinsic curvature and 

(■■■' = e ^k{Ri--M) 
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the Cotton tensor, and the covariant derivatives are de- 
fined with respect to the spatial metric gij. e yfc is the 
totally antisymmetric unit tensor, A is a dimensionless 
constant and the variables k, w and p are constants with 
mass dimensions —1, and 1, respectively. Finally, we 
mention that in action §3§ we have performed the usual 
analytic continuation of the parameters p and w of the 
original version of Hofava-Lifshitz gravity, since such a 
procedure is re quir ed in order to obtain a realistic cos- 
mology [HI, HE IH, Hl| (although it could fatally affect 
the gravitational theory itself). Therefore, in the present 
work A is a positive constant, which as usual is related 
to the cosmological constant in the IR limit. 

In order to add the matter component (including both 
dark and baryonic matter) in the theory one can follow 
two equivalent approaches. The first is to introduce a 
scalar field [29|, |30j and thus attribute to dark matter 
a dynamical behavior, with its energy density p m and 
pressure p m defined through the field kinetic and poten- 
tial energy. Although such an approach is theoretically 
robust, it is not suitable from the phcnomenological point 
of view since it requires special matter-potentials in order 
to acquire an almost constant matter equation-of-state 
parameter (w rn — p m /Pm) as it is suggested by obser- 
vations. In the second approach one adds a cosmolog- 
ical stress-energy tensor to the gravitational field equa- 
tions, by demanding to recover the usual general relativ- 
ity formulation in the low-energy limit [131 . [ill ]. Thus, 
this matter-tensor is a hydrodynamical approximation 
with p m and p m (or p m and w rn ) as parameters. Sim- 
ilarly, one can additionally include the standard-model- 
radiation component (corresponding to photons and neu- 
trinos), with the additional parameters p r and p r (or p r 



and w r ). Such an approach, although not fundamental, 
is better for a phenomenological analysis, such the one 
performed in this work. 

Now, in order to focus on cosmological frameworks, we 
impose the so called projectability condition [ll| and use 
an FRW metric, 



N = l, g t3 = a 2 (t)^ 
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dx l dx J = 
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where K <,=,> corresponding to open, flat, and 
closed universe respectively (we have adopted the con- 
vention of taking the scale factor a(t) to be dimensionless 
and the curvature constant K to have mass dimension 2). 
By varying N and , we obtain the equations of motion: 
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where we have defined the Hubble parameter as H = 
^. As usual, p m (dark plus baryonic matter) follows the 
standard evolution equation 

p m + 3H(p m + Pm ) = 0, (10) 

while p r (standard-model radiation) follows 

p r + 3H(p r +p r ) = 0. (11) 

Finally, concerning the dark-energy sector we can de- 
fine 
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The term proportional to a~ A is the usual "dark radia- 
tion term" , present in Hofava-Lifshitz cosmology [1^, [3(| , 
while the constant term is just the explicit cosmological 
constant. Therefore, in expressions (fTS]> . ([ 13 [) we have 
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defined the energy density and pressure for the effec- 
tive dark energy, which incorporates the aforementioned 
contributions. Finally, note that using (|12|1 . (p~3 j) it is 
straightforward to see that these dark energy quantities 
satisfy the standard evolution equation: 



Pde + 3H(pde + Pde) = 0. 



(14) 



Using the above definitions, we can re- write the Fried- 
mann equations ([8]), ([9]) in the standard form: 



H 
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Therefore, if we require these expressions to coincide with 
the standard Friedmann equations, in units where c = 1 
we set El |33 : 



G 
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(17) 



with G the usual Newton's constant. Note that the run- 
ning of the light speed with A, is not a problem, since 
in this work we will remain in the phenomenologically 
relevant case of A = 1. 



B. Beyond Detailed Balance 

The aforementioned formulation of Hofava-Lifshitz 
cosmology has been performed under the imposition of 
the detailed-balance condition. However, in the litera- 
ture there is a discussion whether this condition leads to 
reliable results or if it is able to reveal the full informa- 
tion of Hofava-Lifshitz gravity (2I [30j] . Therefore, one 
needs to investigate also the Friedmann equations in the 
case where detailed balance is relaxed. In such a case one 
can in general write [U El EE El, SI : 
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where uq = k 2 /12, and the constants <7j are arbitrary 
(with cr 2 being negative) . Note that one could absorb the 
factor of 6 in redefined parameters, but we prefer to keep 
it in order to coincide with the notation of [H, [4l[ . As 
we observe, the effect of the detailed-balance relaxation 
is the decoupling of the coefficients, together with the 
appearance of a term proportional to a~ 6 . In this case the 
corresponding quantities for dark energy are generalized 
to 

= 21 j- azK2 

~ 6 + 
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3 ~ 6 + 18a 4 



PDE\ 



Pde\ 



6a 4 
^K 2 



a 4 K 



(20) 
(21) 



Therefore, it is easy to see that 

PcsU-db + ZH{ PDE \ aon _ dh +p DE L n . dh ) = 0. (22) 

Finally, if we force (fT8|) , (fT9|) to coincide with the standard 
Friedmann equations, we result to: 

6cr 
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III. OBSERVATIONAL CONSTRAINTS 

Having presented the cosmological equations of a uni- 
verse governed by Hofava-Lifshitz gravity, both with and 
without the detailed-balance condition, we now proceed 
to study the observational constraints on the model pa- 
rameters. This is performed in the following two subsec- 
tions, for the detailed and non-detailed balance scenarios 
separately. We mention that since the cosmological ob- 
servations lie deep inside the IR, in the following we set 
the running parameter A to 1. 



A. Constraints on Detailed-Balance scenario 

We work in the usual units suitable for observational 
comparisons, namely setting %ixG — 1 (we have already 
set c = 1 in order to obtain ([17])). This allows us to 
reduce the parameter space, since in this case (fiTf lead 
to: 



p 2 A = 2. 

Inserting these relations into Friedmann equation 
obtain 
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In terms of the usual density parameters (O m = 
p m /(3H 2 ), n K = -K/(H 2 a 2 ), n r = p r /(3H 2 )) this ex- 
pression becomes: 
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Finally, applying this relation at present time and setting 
the current scale factor clq — 1 we get: 

1 (K 2 A\ 
1 - Q m Q — fl r0 - Qkq = 2A + ~2 ) ' 

where a O-subscript denotes the present value of the cor- 
responding quantity. 

As was mentioned above, we have used the analytic 
continuation, as a result of which A is positive. Thus, 
relation (|2"T|) can in principle be satisfied by a suitable 
choice of A. However, note that without the analytic 
continuation (and therefore with a negative A) relation 
(|2"7[) could never be satisfied (as expected, since in this 
case the theory would not have the A = 1 IR, general- 
relativity limit), and this offers another indication, from 
the phenomenological point of view, for the necessity of 
the analytic continuation in the detailed-balance version 
of Hofava-Lifshitz cosmology. 

In order to proceed to the elaboration of observa- 
tional data, we consider as usual the matter (dark plus 
baryonic) component to be dust, that is w m « 0, and 
similarly for the standard-model radiation we consider 
w r = 1/3, where both assumptions are valid in the epochs 
in which observations focus. Therefore, the correspond- 
ing evolution equations (fTOj) . (|1 1[) give p m = p m o/a 3 and 
Pr = Pro/° 4 respectively. Finally, instead of the scale 
factor it proves convenient to use the redshift z as the in- 
dependent variable, which is given by X+z = a^/a = I /a. 
Inserting these into Friedmann equation (|25p we obtain 



H 2 



Hi [n m0 {\ + zf + n r0 (i + zf + n K0 (i + zf + 

L (l + ^) 4 ]}, (28) 
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where we have also introduced the dimcnsionlcss param- 
eter u = A/(2H 2 ). Thus, the constraint (f27|) can be 
rewritten as: 



n 
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As we have already mentioned above, the term 
n 2 K0 /(Auj) is the coefficient of the dark radiation term, 
which is a characteristic feature of the Hofava-Lifshitz 
gravitational background. Since this dark radiation com- 
ponent has been present also during the time of nucle- 
osynthesis, it is subject to bounds from Big Bang Nucle- 
osynthesis (BBN). As discussed in more details in the 
Appendix, if the upper limit on the total amount of 
dark radiation allowed during BBN is expressed through 
the parameter AN ^ of the effective neutrino species 
PM [lOi [lOl [Toll, then we obtain the following con- 
straint : 



o2 

— ^ = 0.135AA^fi r0 . 
4w 



(30) 



Finally, we mention that as usual, the density parameter 
for standard model radiation (photons and three species 



of neutrinos) £! r o is entirely determined by f^mf), Hq and 
the measured value of the CMB temperature [83] . 

In most studies of dark energy models it is customary 
to ignore curvature (e.g.JH M, M,M,M,M,M,^), 
especially concerning observational constraints. This 
practice is well motivated for at least two reasons. 
Firstly, most inflationary scenarios predict a high de- 
gree of spatial flatness. Secondly, the CMB data im- 
pose stringent constraints on spatial flatness in the con- 
text of constant-it; models (for example a combination 
of WMAP+BAO+SNIa data [87[ provides the tight si- 
multaneous constraints —0.0179 < £Iko _ 0.0081 and 
-0.12 < 1 + w < 0.14, both at 95% confidence). 

However, it is important to keep in mind that ow- 
ing to degeneracies in the CMB power spectrum (see 
[961 ] and references therein), the limits on curvature de- 
pend on assumptions regarding the underlying dark en- 
ergy scenario. For example, if instead of a constant w 
one works with a linearly varying w, parameterized as 
w (a) — wo + (1 — a) w a , the error on Q r-q is much larger, 
on the order of a few percent [97|, [H, [9!| . The constraints 
on curvature for different parameterizations was studied 
in [lQ^, [lOll, Eol] . The authors of [lO^ showed that for 
some models of dark energy the constraint on the curva- 
ture is at the level of 5% around a flat universe, whereas 
for others the data are consistent with an open universe 
with VLko ~ 0.2. According to [§8|, geometrical tests 
such as the combination of the Hubble parameter H(z) 
and the angular diameter distance Da(z), using (future) 
data up to sufficiently high redshifts z ~ 2, might be 
able to disentangle curvature from dark ener gy eyolu- 
tion, though not in a model-independent way. [1031 1 104| 
highlight the pitfalls arising from ignoring curvature in 
studies of dynamical dark energy, and recommend treat- 
ing Qko as a free parameter to be fitted along with the 
other model parameters. 

In the present work, the spatial curvature plays a very 
crucial role, since, as it has been extensively stated in 
the literature [29|, [§3], Hofava-Lifshitz cosmology coin- 
cides completely with ACDM if one ignores curvature. 
Therefore, and following the discussion above, we choose 
to treat Qko as a free parameter. 

In summary, the scenario at hand involves four pa- 
rameters (we fix Hp by its 5-year WMAP best-fit values, 
given in Table 1 of [83]), namely fl m o, Qko, w and A.N V , 
subject to constraint equations ([29)) and ([30| . Therefore, 
only two of these parameters are independent. Although 
one would usually expect to be able to choose two of them 
at will, the non-linear nature of the constraint equations 
does not facilitate this, and one has no choice but to use 
f2 m o and AA?"„ as free parameters. Inverting ([2"9"[) and 
(|30p to express ui and flxo m terms of these independent 
parameters for a given curvature, we obtain: 

lu (K; n ro0 , AN U ) = 1 - Q m0 - (1 - A^)^ 



0.73 sgn (k) VAiv, J n r0 - n mQ n r0 - n 2 r0 (31) 
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and 



\n K0 (O m0 , AN V ) | = v/0.54 AJV„ n r0 w (fi m o, AJV„). 

(32) 

The BBN upper limit on AiV^ is -1.7 < AiV„ < 2.0 
[IM EM EM Ell- A negative value of AN V (which is 
usually associated with models involving decay of a mas- 
sive particle) is not possible in the present model, since uj 
(i.e. A) is always positive. Additionally AN V = corre- 
sponds to the zero curvature scenario (a non-interesting 
case since Hof ava-Lifshitz cosmology with zero curvature 
becomes indistinguishable from ACDM). 

In Fig. Q] and Fig [5] we use a combination of obser- 
vational data from SNIa, BAO and CMB to construct 
likelihood contours for the parameters £l m o and AN V for 
positive and negative curvatures respectively. These fig- 
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FIG. 1: (Color Online) Contour plots of Q,m0 vs AJV„ for 
positive curvature (K > 0), under SNIa, BAO and CMB ob- 
servational data. The yellow (light) region is excluded at the 
2a level, and the orange (darker) region is excluded at the la 
level. The red (darkest) region is not excluded at either confi- 
dence level. The white diamond marks the best-fit point. The 
model parameters u = A/(2Hq) and Qko = ~K/{Hq) are 
related to f2 m o and AN V through $Tj) and 
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FIG. 2: (Color Online) Contour plots of Q m o vs AN V for 
negative curvature (K < 0), under SNIa, BAO and CMB 
observational data. The yellow (light) region is excluded at 
the 2a level, and the orange (darker) region is excluded at the 
la level. The red region is not excluded at either confidence 
level. The white diamond marks the best-fit point. The model 
parameters uj = A/(2Hq) and Qko = —K/(Hq) are related to 
fl m o and AN V through (3l\) and KM) . 



In conclusion, we have shown that Hofava-Lifshitz cos- 
mology under the assumption of detailed-balance con- 
dition cannot fulfill observational requirements without 
the analytic continuation transformation. Under analytic 
continuation the observational constraints on the param- 
eters are quite tight. This feature was already mentioned 
in [l5j |. following qualitative theoretical arguments con- 
cerning the effective light speed in Hofava-Lifshitz frame- 
work, where it was stated that a fine tuning would be 
needed as a way out. The analysis of this section offers 
new, phenomenological indications towards the direction 
of tight constraints. 

B. Constraints on Beyond-Detailed-Balance 
scenario 



ures show that the Hofava-Lifshitz cosmological scenario 
under the detailed balance condition is not ruled out by 
observations. However, they lead to tight constraints on 
the amount of dark radiation allowed at the time of nu- 
cleosysnthesis (tighter than the corresponding limits from 
BBN), and thus to the parameter A. For example, the 
la limits on £Iko, A and fi (which is actually connected 
to A through (|17p) are presented in Table HI For conve- 
nience we have kept the factors of Hq and 8ttG. Thus, 
one can either use the usual ansatz 8irG = Hq = 1, or 
insert physical units using H = 1.503 x 10~ 42 GeV and 
8ttG = 1.681 x 10~ 37 GeV" 2 . In the later case, one ob- 
tains < A < 1.86 x 10" 9 eV 4 and 3.30 x 10 60 < \i < oo 
for the positive curvature case, and similarly for the neg- 
ative curvature one. 



In units where 8irG = 1 relations 

tfo = 1/3 
o- 2 = -6. 



give 



(33) 



Using these values and following the procedure of the 
previous subsection, the Friedmann equation (|18|) can be 
written as 



H 1 = H, 



jjOno (i + zf + n r0 (i + zf + n K0 (i + zf + 

+ ui + u; 3 (1 + zf + lu 4 (1 + zf }. (34) 

In this expression we have introduced the dimensionless 
parameters u>±, lo% and w±, related to the model parame- 
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{9,-kG/HI) A 




4 


(0, 0.0038) 


(0, 1.4189) 


(1.1872, oo) 


4 


(-0.0039, 0) 


(0, 1.4063) 


(1.1925, oo) 



TABLE I: la limits on the parameter values for the detailed-balance scenario, for positive and negative curvature. 



ters <Ti, <T3 and (74 through: 



CJ4 
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<t 3 H%S1 2 K0 
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(35) 



Furthermore, we consider the combination 074 to be posi- 
tive, in order to ensure that the Hubble parameter is real 
for all redshifts. uj 4 > is required also for the stability of 
the gravitational perturbations of the theory For 
convenience we moreover assume 03 > 0, that is lu 3 > 0. 

The scenario at hand involves the following parame- 
ters: Hq, Q m o, fluo, wi, L03 and 0J4. Similarly to the 
detailed-balance section these are subject to two con- 
straints. The first one arises from the Friedman equation 
at z = 0, which leads to 



O m0 + f2 r o + &K0 + Ml + W3 + OJ4, = 1. 



(36) 



This constraint eliminates the parameter w\. 

The second constraint arises from BBN considera- 
tions. The term involving lu 3 represents the usual dark- 
radiation component. In addition, the W4-term represents 
a kination-like compone nt (a qui ntessence field domi- 
nated by kinetic energy |l09l . IllOl ]). If AN V represents 
the BBN upper limit on the total energy density of the 
universe beyond standard model constituents, then as we 
show in the Appendix we acquire the following constraint 



at the time of BBN (z = z B bn) [Ml GM M 



g const r. 
a [III: 



lj 3 + uj a (1 + z B bn) = 0.135AA^ 



HI- 



(37) 



It is clear that BBN imposes an extremely strong con- 
straint on 0J4, since its largest possible value (correspond- 
ing to lu 3 = 0) is ~ 10" 24 . 

We use relation (|37[) to eliminate 104 in favor of L03 and 
AN Ul and treat lu 3 and £Iko a s our free parameters. Since 
(x>4 is non-negative, relation (|37p determines also the up- 
per bound of u>3. For the remaining parameters, £l m o 
and Hq (unless otherwise stated) we assume priors cor- 
responding to their 5-year WMAP best-fit values (given 
in Table 1 of [13). 

We now proceed to constrain the free parameters 
Qko and LO3 through observations. Using the combined 
SNIa+CMB+BAO data, we construct likelihood con- 
tours for these two parameters. The details and the tech- 
niques of the construction are presented in the Appendix. 
Furthermore, since the BBN limits on AN V vary depend- 
ing on assumptions, in addition to our canonical choice 
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FIG. 3: (Color Online) Contour plots of log 10 (u>3) vs 
log 10 |ftifo| for K > and AN U = 0.1, using SNIa, BAO 
and CMB data. The yellow (light) region is excluded at the 
2a level, and the orange (darker) region is excluded at the 
la level. The red (darkest) region is not excluded at either 
confidence level. The white diamond marks the best-fit point. 
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FIG. 4: (Color Online) Contour plots of log 10 (u>3) vs 
log 10 |fijro| for K < and AN„ = 0.1, using SNIa, BAO 
and CMB data. The yellow (light) region is excluded at the 
2a level, and the orange (darker) region is excluded at the la 
level. The red (darkest) region is not excluded at either confi- 
dence level. The white diamond marks the best-fit point (near 
the top right corner in this case). 
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a /(8nG) 






(8-KG/Hi) a x 


(72 




(7 4 /(8ttG) 


1/3 


0.1 


(0, 0.01) 


(4.29, 4.33) 


-6 


(0, 0.03) 


(-9.08 x 10~ 22 , 0) 


1/3 


0.1 


(-0.01, 0) 


(4.40, 4.45) 


-6 


(0, 0.81) 


(0, 5.66 x 10~ 22 ) 


1/3 


2.0 


(0, 0.04) 


(4.13, 4.45) 


-6 


(0, 0.01) 


(-1.77 x 10~ 20 , -2.62 x 10 " 21 ) 


1/3 


2.0 


(-0.01, 0) 


(4.40, 4.45) 


-6 


(0, 0.23) 


(-2.61 x 10~ 20 , -1.16 x 10" 20 ) 



TABLE II: lcr limits on the parameter values for the beyond-detailed-balance scenario, for positive and negative curvature, and 
for two values of the effective neutrino species parameter AJV„ (see text). 




iog 10 (w3) 

FIG. 5: (Color Online) Contour plots of log 10 (u>3) vs 
log 10 |fiA-o! f°r K > and AN„ = 2.0, using SNIa, BAO 
and CMB data. The yellow (light) region is excluded at the 
2a level, and the orange (darker) region is excluded at the 
la level. The red (darkest) region is not excluded at either 
confidence level. The white diamond marks the best-fit point. 




FIG. 6: (Color Online) Contour plots of log 10 (u>3) vs 
log 10 \tt K -o\ for K < and AN„ = 2.0, using SNIa, BAO 
and CMB data. The yellow (light) region is excluded at the 
2a level, and the orange (darker) region is excluded at the 
la level. The red (darkest) region is not excluded at either 
confidence level. The white diamond marks the best-fit point. 



of upper limit AA^ = 2.0, we have also considered the 
more restrictive limit of AN V = 0.1. 

Figures [3] and |4] depict the la and 2a U3 — \£Iko\ con- 
tours, for AiV„ = 0.1, for positive and negative curvature 
respectively. Figures [5] and [5] are the corresponding plots 
using AN V = 2. In each case, CJ3 extends over its entire 
allowed range, namely < 0J3 < 0.135AN v n r Q. 

As we observe, the Hofava-Lifshitz cosmological sce- 
nario beyond the detailed balance condition is not ruled 
out by observations. However, they impose strong con- 
straints on W3 (for the case of AN V — 2.0 the constraints 
on UJ3 for both the positive and negative curvature cases 
are stronger than the upper bound from BBN), and ex- 
tremely tight constraints on LJ4. Additionally, the con- 
straints on the curvature are of the order of a percent. 
Note that the contours expand as AN V is reduced. This 
is expected since the smaller the amount of exotic compo- 
nents (dark-radiation and kination-like ones), the closer 
the model is to ACDM. 

The approximate la limits on the model parameters 
<7j are presented in Table [TTJ The lower limit on a 3 is 
zero. From (|3"6")l and (|37p it is clear that a\ and (74 at- 



tain their maximum values when (73 is at its minimum, 
and vice versa. Similarly to the previous subsection, one 
can either use the usual ansatz 8irG = Hq = 1, or in- 
sert physical units using H Q — 1.503 x 10~ 42 GeV and 
8ttG = 1.681 x 10~ 37 GeV" 2 . In the later case, one 
obtains that a\ is tightly constrained to be at the level 
of the cosmological constant (10~ 12 eV 4 ), as expected. 
Additionally, the data impose extremely stringent con- 
straints on (74, which was also expected. However, even 
for such small values, the phenomenological implications 
of the kination-like (74-co mponent a re very interesting. 
As discussed in detail in jllll . Ill2j . it could dominate 
the universe prior to BBN and it could significantly af- 
fect the freeze-out, and hence the relic abundances of 
neutralino dark matter, by a few orders of magn itude . 
For dark matter that decays into leptons (see e.g. jll3l |) 
this could be relevant to recent observations o f high en- 
ergy positr ons a nd electrons by the PAMELA [l!4l 1 1 1 5| 
and ATIC [116| | experiments. 
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IV. CONCLUSIONS 

In this work we constrained Hof ava-Lifshitz cosmology 
using observational data. In particular, we considered 
scenarios where the gravitational sector is forced to sat- 
isfy the detailed-balance condition, and also those where 
this condition is relaxed. Additionally, we have included 
the matter and radiation sectors following the usual ef- 
fective fluid approach. These constructions, which cover 
the range of Hof ava-Lifshitz cosmology, were confronted 
with data from BBN, SNIa, CMB and BAO observations. 

Our first result is that the detailed-balance formu- 
lation of Hofava-Lifshitz gravity cannot fulfill observa- 
tional requirements, without the analytic-continuation 
transformation. Under the analytic continuation we 
found that Hofava-Lifshitz cosmology can be compatible 
with observations, and we presented the corresponding 
contour-plots on the model parameters. These likelihood- 
contours impose tight constraints on the model param- 
eters, and the corresponding lcr-bounds are presented 
in Table Q] However, we mention that although ana- 
lytic continuation is necessary for a realistic cosmology, it 
can fatally affect the gravitational theory itself, spoiling 
its initial stability and well-behaving nature. Therefore, 
the detailed-balance version of Hofava-Lifshitz cosmol- 
ogy seems rather unlikely to be a robust description of 
nature. 

The version of Hofava-Lifshitz cosmology in which 
the detailed-balance condition has been abandoned, is 
also compatible with observations. We constructed the 
likelihood-contours for the two involved free parame- 
ters, namely the curvature and the dark-radiation co- 
efficients. As we showed, observations lead to strong 
bounds in these parameters, and the corresponding la- 
allowed ranges are presented in Table HH This feature 
was expected, since the data refer to redshifts in which 
the novel terms of Hofava-Lifshitz cosmology are down- 
graded. However, these terms can have significant cosmo- 
logical implications prior to nucleosynthesis, which could 
be probed by recent observatio ns of high energy posi trons 
and electrons by the PAMELA [lli,lll5j and ATIC flify 
experiments. 

Although the present analysis indicates that Hofava- 
Lifshitz cosmology can be compatible with observa- 
tions, it does not enlighten the discussion about possible 
conceptual problems and instabilities of Hofava-Lifshitz 
gravity, nor it can interfere with the questions concern- 
ing the validity of its theoretical background, which is the 
subject of interest of other studies. In particular, without 
a solid theoretical basis, it is not clear whether Hofava- 
Lifshitz gravity is able to pass the basic parametrized 
post newtonian (PPN) tests that any physically inter- 
esting gravitational theory should jll7l . Ill8l lll9L 1 120j ] . 
The present work just faces the problem from the phe- 
nomenological point of view, and thus its results can been 
taken into account only if Hofava-Lifshitz gravity passes 
successfully the aforementioned theoretical tests. 
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APPENDIX: OBSERVATIONAL DATA AND 
CONSTRAINTS 

In this appendix, we briefly review the main sources 
of observational constraints used in this work, namely, 
Big Bang Nucleosynthesis (BBN), Baryon Acoustic Os- 
cillations (BAO) and the Cosmic Microwave Background 
(CMB). 

a. Big Bang Nucleosynthesis constraints 

Big Bang Nucleosynthesis (BBN) provides a highly 
sensitive tool for probing p hysi c s beyond the standard 
model (for reviews see e.g. [lOl DM DM [Tol] ) . Abun- 
dances of light elements predicted by BBN, particularly 
the 4 He one, are sensitive to the expansion rate of the 
universe (or equivalently to its total energy density) at 
the time of BBN. Additionally, the abundances depend 
also on the baryon to photon ratio, though this ratio can 
be independently determined from the CMB by WMAP 
data [87| . BBN therefore imposes constraints on the den- 
sities of possible extra exotic radiation constituents (be- 
yond the standard model photons and three flavors of 
neutrinos). 

The constraints on the energy density of these exotic 
constituents are usually expressed in terms of the effec- 
tive neutrino species AN U . Assuming that neutrinos are 
fully decoupled from photons and do not gain energy from 
annihilation, they are colder than photons by a factor 
of T„/T 7 = (4/lf) 1/3 (see e.g. [l07|). Using also that the 
neutrino energy is related to the photon one by a factor 
of 7/8, the total energy density of relativistic species - 
photons and (3 + A A'',,) species of neutrinos - reads 

PT = Pl + (3 + AN U ) f i.) ' p y . (A.l) 

In terms of the total standard-model relativistic density 
(photons and three species of neutrinos) /? 7 „, the above 
can be written as 

p T = (1 + 0.135AAy p lv . (A.2) 

In the pre sent work we use the upper limits on AN V 
provided in |l07| : -1.7 < AN V < 2.0, although more 
restr i ctive bounds have also been proposed (see e.g. 
[l2ll [l22]). We mention that the above limits do not 
apply to models in which the dark radiation or other 
exot ic components are injected later than BBN (see 
[l23| for an example of such a model), however they are 
obviously applicable to Hofava-Lifshitz cosmology, in 
which dark radiation is always present, arising from the 
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gravitational theory itself. 

b. Type la Supernovae constraints 

In order to incorporate supernova co nstra ints we use 
the Union08 compilation of Snla data |l24j . This is a 
heterogeneous data-set, consisting of data from the Su- 
pernova Legacy Survey, the Essence survey, the recently 
extended data-set of distant supernovae observed by the 
Hubble Space Telescope, as well as older data-sets. 

The x 2 from SNIa is calculated as: 



N 

J2 [Mobs (zi) ~ thh {zi)Y 

2 _ i=l 



(A.3) 



where N = 307 is the number of SNIa data points. ^ b s is 
the observed distance modulus, defined as the difference 
between the apparent and absolute magnitude of the su- 
pernova. The a^i are the errors in the observed distance 
moduli, arising from a variety of sources, and assumed to 
be gaussian and uncorrelated. The theoretical distance 
modulus fith depends on the model parameters a, via the 
dimensionless luminosity distance Di(z;ai): 



PZ TT 

D L (z; Oi) = (1 + z) / dz' 

Jo H(z';a,i) 



(A.4) 



as follows: 



/i th (z) = 42.38 - 5 log 10 h + 5 log 10 [D L (z; en)] . (A.5) 

The marginalization over the present value of the Hubble 
parameter is p erformed following the techniques de- 
scribed in 125], and we construct x 2 likelihood contours 
for the various model parameters. 

c. CMB constraints 

We use the CMB data to impose constraints on the 
parameter space, following the re cipe described in [87j ]. 
The "CMB shift parameters" [si, [l26| are defined as: 

R = \/n m0 H r (z*) , l a = Trr (z») jr s (z*) . (A. 6) 

R can be physically interpreted as a scaled distance to 
recombination, and l a can be interpreted as the angular 
scale of the sound horizon at recombination. r{z) is the 
comoving distance to redshift z defined as 



r(z) 



1 



o H(z) 



dz, 



(A.7) 



while r s (z„) is the comoving sound horizon at decoupling 
(redshift z*), given by 



, H(z)y/3(l + Rb/(l + z)) 



zdz. 



(A.8) 



The quantity R^ is the ratio of the energy density of pho- 
tons to baryons, and its value can be calculated as Rb = 



31500ft b0 /i 2 (T C ub/2.7K) 4 , (O 60 being the present day 
density parameter for baryons) using Tcmb = 2.725 [871 ]. 
The redshift at decoupling (Qbo, ii m n, h) c an be calcu- 
lated from the following fitting formula |l27j : 



1048 



0.00124 (n b0 h 2 ) 



-0.738 



1 + gi (fl m oh 2 ) 



2\92 



with gl and g2 given by: 



(A.9) 



-0.238 



.9i = 



92 = 



0.0783 (n bQ h 2 ) 

1 + 39.5 (tt b0 /i 2 ) - 763 
0.560 



i + 2i.i(n M ^ 2 ) 1 " 81 ' 

Finally, the \ 2 contribution of the CMB reads 

XCMB = VcMB C invV CM B- (A. 10) 

Here V C mb = P-Pdata, where P is the vector (l a ,R, z*) 
and the vector Pdata is formed from the WMAP 5-year 
maximum likelihood values of these quantities [87] ]. The 
inverse covariance matrix Cj nv is also provided in [87j ■ 

d. Baryon Acoustic Oscillation constraints 

In this case the measured quantity is the ratio d z = 
r s (zd) /Dy (z), where Dy (z) is the so called "volume 
distance" , defined in terms of the angular diameter dis- 
tance Da = r (z) I (1 + z) as 



D v (z) 



{l + zfD\ {z)z 
H(z) 



-, 1/3 



(A.ll) 



and Zd is the redshift of the baryon drag epoc h, which 
can be calculated from the fitting formula [1281 ]: 



1291 (n m0 h 2 



0.251 



1 + (fW 2 )°- 828 

where b\ and hi are given by 



h (n b0 h 2 



(A.12) 



h = 0.313 (n m0 h 2 ) °' 419 1 + 0.607 (n m0 h 2 ) 



2x0.674 



62 = 0.238 (nmoh" 



,0.223 



We use the two m easurements of d z at redshifts z = 0.2 
and z = 0.35 [l29j . We calculate the \ 2 contribution of 
the BAO measurements as: 



XBAO — VB AO C inv V B AO- 



(A.13) 



Here the vector Vbao = P Pdata, with P = 
(do.2,do.3s), and P da ta = (0.1905, 0.1097), the two mea- 
sured BAO data poi nts [ l29t ] . The inverse covariance 
matrix is provided in [1291 ] . 
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